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BOUNDARY CONDITIONS 



The colliding plane wave problem has been formulated in the previous 
chapters as a characteristic initial value problem with initial data specified 
on two null hypersur faces. At this point it is necessary to discuss the 
junction conditions between the four regions described in Figure 3.1, and 
particularly on the initial boundaries of region IV. 

7.1 General discussion 

When pasting together particular solutions of Einstein's equations ap- 
plying to different regions of space-time, it has generally been thought 
that the appropriate junction conditions are those of Lichnerowicz (1955). 
These require that there exist coordinates in which the metric tensor is 
C 1 and piecewise at least C 2 . These conditions appear to be reasonable 
since the curvature tensor involves the second derivatives of the metric 
tensor. 

Work on junction conditions up to 1966 has been surveyed by Israel 
(1966), with particular emphasis on conditions across non-null bound- 
aries. He has also shown that, if the distributional part of the energy- 
momentum tensor is zero, then the second fundamental forms induced by 
the metric match across the boundary, and this implies that the full Rie- 
mann tensor is regular. This condition was proposed by Darmois (1927), 
and has been shown to be equivalent to the Lichnerowicz conditions by 
Bonnor and Vickers (1981). 

An alternative set of conditions has been proposed by O'Brien and 
Synge (1952). In the non-null case, these also have been shown by Israel 
(1958) and Robson (1972) to be equivalent to those of Lichnerowicz. How- 
ever, since these are not stated covariantly, some differences may occur 
as discussed by Bonnor and Vickers (1981). 

O'Brien and Synge (1952) have also proposed conditions across null 
boundaries that are weaker than those of Lichnerowicz. Denoting the null 
hypersurface by x° = constant, with goo = 0, they propose simply that 
the components 

9i*u, 9 tJ 9ij,o, 9 l °9ij,o (7.1) 
where = 1,2,3), should be continuous across x° = constant. 
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By requiring that the metric tensor be C and piecewise at least 
C 2 , the Lichnerowicz conditions specifically exclude the impulsive grav- 
itational waves that are discussed in Chapter 3. It can be shown, how- 
ever, that the line element (3.3) does satisfy the O'Brien-Synge conditions 
across the null hypersurface u = 0. This connects the metrics (3.6) and 
(3.7) in regions I and II. The importance of an appropriate choice of co- 
ordinates may also be pointed out here, since, when the continuous line 
element (3.3) is transformed to the form (3.1), the metric tensor then 
contains a delta function discontinuity. 

It has been shown by Robson (1973), that the appropriate junction 
conditions across a null hypersurface are those of O'Brien and Synge 
(1952). Junction conditions across null hypersurfaces have been further 
analysed by Penrose (1972) and Clarke and Dray (1987). In addition, Bell 
and Szekeres (1974) have shown that, for colliding plane electromagnetic 
waves, the Lichnerowicz conditions have to be relaxed in favour of those 
of O'Brien and Synge. 

7.2 Junction conditions for colliding plane waves 

In order to discuss the collision of plane waves, it has been found conve- 
nient to divide space-time up into four regions as described in Figure 3.1. 
These regions are bounded by the two null hypersurfaces u = and v = 0. 
It has also been suggested in Section 6.4, that it is convenient to assume 
that the line element (6.20) applies to the entire space-time. However, 
the metric functions U, V, W and M must take different forms in the 
four regions. 

It is now necessary to consider the conditions that should be imposed 
on these functions at the boundaries of the four regions. The O'Brien- 
Synge conditions (7.1) imply that V, W and M are continuous and that 
U is smooth across these null boundaries. However, according to (6.24), 



It is therefore necessary that f(u) and g(v) are at least C 1 . 

In order for the line element (6.20) to describe a collision of plane 
waves, it must be assumed that U, V, W and M are functions of u and v 
in region IV, are functions of u only in region II, are functions of v only 
in region III, and are constants in region I. Concentrating initially on the 
metric function e~ u = f + g, it is appropriate to choose 



U = -\og(f(u) + g(v)). 



(7.2) 
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for u < 0, 
for v < 0, 



f (0) = 
g'(0) = 0. 



(7.3) 
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The approaching waves are then partly described by f(u) in region II, 
and by g(v) in region III. With this choice it may now be noted that 

in region I : e~ u = 1 

in region II: e~ u = \ + f(u) 

(7.4) 

in region III : e~ u = \ + g(v) 
in region IV : e~ u = f(u) + g(v). 

The junction condition that U is smooth across the boundaries requires 
that the function f(u) must have the same form in both regions II and 
IV. Similarly g(v) must have the same form in regions III and IV. 

At this point it is convenient to concentrate on the initial boundaries 
between regions I and II, and between regions I and III. It is then possible 
to use the transformations (6.7) to put M = in all three of these initial 
regions. Equations (6.226, c), which now apply only to regions III and II 
respectively, become 



g" = \e u g' 2 - \e~ v \V V 2 cosh 2 W + W 2 + 4$°$°) 
/" = \e u f 2 - \e~ u (y 2 cosh 2 W + W 2 U + 4d>°$°). 



(7.5) 



These equations, together with (7.3), imply that / and g are monoton- 
ically decreasing functions for positive arguments, at least in regions II 
and III. The scale transformations (6.7) cannot reverse the directions of 
the parameters u and v. Thus, even in region IV with M non-zero, / and 
g must be monotonically decreasing functions. 

It is always possible, therefore, to use the transformation (6.7) to 
express the functions / and g in the forms 

/ = i _ ( ClU )»ie(«), g = \- {c 2 v) n *Q{v) (7.6) 

which apply globally. It is also possible to use further transformations to 
put ci = 1 and c 2 = 1. However, these particular constants describe a 
measure of the magnitudes of the approaching waves, and it is therefore 
often convenient to retain them without this further rescaling. For situa- 
tions in which M = in the inital regions I, II and III, the product c\c 2 
has an invariant meaning as the amount of non-linearity in the interac- 
tion. However, if transformations are used to put c\ = 1 and c 2 = 1, then 
this magnitude is absorbed into the metric component e~ M . 

The expressions (7.6) may be of particular convenience, as they trans- 
parently demonstrate the technique by which an exact solution obtained 
in region IV in terms of these functions can be extrapolated back to re- 
gions III, II and I. However, this particular parametrization is not always 
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the most convenient, as will be demonstrated for example in the Bell- 
Szekeres solution to be described in Chapter 15. 

It may also be noted that in the line element (6.20) the metric coeffi- 
cient e _c/ , according to (6.24), is given by f(u) +g(v) where / and g have 
now both been shown to be decreasing functions from the value |. It is 
therefore inevitable that a singularity will develop as / + g — > 0. Whether 
this is a curvature singularity or simply a coordinate singularity will have 
to be considered in detail. An initial discussion of this topic is given in 
the next chapter. 

In practice, it is not easy to find exact solutions in the interaction 
region for any specified initial conditions in which the metric functions in 
regions I, II and III are the initial data. Although this approach will be 
discussed later, most of the explicit exact solutions that will be described 
in the following chapters have been obtained by first solving the field 
equations in region IV, and then extrapolating back to determine the 
approaching waves that would give rise to them. In this context it may 
be noted that, since / and g are monotonically decreasing functions, they 
may be adopted as coordinates in the interaction region. Also, since these 
functions are required to be smooth across the boundaries, they may easily 
be matched to the null coordinates u and v in regions II and III. 

For colliding gravitational waves, attention is concentrated on equa- 
tions (6.22<i, e). In the interaction region it is possible to use / and g as 
coordinates and these equations may then be integrated to give V(f,g) 
and W{f, g) subject to the initial data given by V{\,g) and W(\,g) on 
the surface u = and by V(/, |) and W(f, \) on the surface v = 0. The 
situation for colliding electromagnetic waves is a little more complicated 
with the addition of (6.21) as well as the extra terms in (6.22d, e). How- 
ever, even in this case the junction conditions for V and W do not usually 
impose additional complications. 

Once V and W are found from equations (6.22<i, e), it is then neces- 
sary to integrate (6.226, c) to obtain M. Because the integrability condi- 
tions are automatically satisfied, such a solution is known to exist. 

Equations (6.226, c) can now be written in the form 



M v = ~j + - cosh 2 W + Wl + 4d>°$°) 

Mu = ~J + 2{JTgj ~ COsh2 W + W "+ **°2*2) ■ 



(7.7) 



It is thus convenient to put 
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where S satisfies 



S g = -k(f + 9) cosh2 W + W l + 



S f = -l(f + g)( Vf cosh 2 W + Wj + 4-^ 



(7.9) 



The boundary conditions discussed above require that e~ be con- 
tinuous, and that /(0) = |, /'(0) = 0, #(0) = | and g'(0) = 0. However, 
in view of (7.8), e~ M cannot be continuous across u = and v = unless 
e~ s is unbounded on these boundaries. For the junction conditions to 
be satisfied, it is therefore essential that the functions V, W, $g an d 3>2> 
which are solutions of equations (6.22<i, e) and (6.21), should be such that 
the solution of (7.9) for S must contain terms of the form 

S = h log(| -f)+k 2 log(i -g) + log(pip 2 ) + • • • (7.10) 

where pi and ^2 are constants. If the leading terms in a power series 
expansion for / and g in region IV have the form 

/=§- (ci«) m + • • • , <7=|- M) n2 + • • • (7-11) 

then e _M is continuous across the boundaries only if S contains the terms 
in (7.10) where the constants k\ and k% are given by 

fci = l-l/m, k 2 = l-l/n 2 . (7.12) 

In order to satisfy (7.3), it is necessary that n\ > 2 and n 2 > 2. It 
therefore follows that fci and fc 2 must be restricted to the range 

\<k u k 2 <l. (7.13) 

In addition, in order to put M zero in region I, it is also appropriate to 
put 

pi = -nici, p 2 = -^2C 2 - (7.14) 

It can be seen that the boundary conditions that are appropriate for 
colliding plane waves become fairly complicated. It is mainly a question 
of choosing solutions V, W ', $q an d ^* 2 of equations (6.22<i, e) and (6.21) 
such that the solution M of (6.226, c) contains terms of the form (7.10). 
Even then, the restrictions (7.12) are placed on the parameters and hence 
on the forms of / and g. 
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Since the main field equations involve the functions V, W, $q and 
$2? it is convenient to consider a form of the boundary conditions that 
applies to these functions only. This can be achieved by substituting 
(7.10) into (7.9) and considering boundaries as / — > 1/2 and g — > 1/2. 
The resulting conditions can be conveniently expressed in the form 



lim 

g^l/2 

lim 

/-1/2 



a-g)lV q 2 cosh 2 W + W"+4 nn 
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\ - /) V/ cosh 2 W + + 4 



/ 



/2 



= 2k 2 
= 2/ci 



(7.15) 



where k\ and fc 2 must satisfy (7.13). It is sometimes more convenient 
to revert to expressions involving the null coordinates u and v. These 
become 



lim 



lim 



V 2 cosh 2 W + Wl + 4%% 

v n 2 —2 

V 2 cosh 2 W + Wl + 4$^ 



u 



ni-2 



= 2n 2 (n 2 - l)c™ 2 
= 2ni(m - l)c^. 



(7.16) 



This form of the conditions is particularly convenient when n\ = n 2 = 2. 

These conditions will be discussed further for colliding gravitational 
waves in Section 11.2 and for the collision of a mixture of electromagnetic 
and gravitational waves in Section 16.1. 

Before concluding this chapter, it is appropriate to consider the possi- 
bilities that arise when the continuity conditions for U are slightly relaxed, 
requiring only that it is continuous so that ni, n 2 > 1. In this case, impul- 
sive components occur in the Ricci tensor on the boundaries between the 
different regions. It is then possible that such components could be inter- 
preted in terms of impulsive null matter fields. Explicit solutions of this 
type have in fact been given by Dray and 't Hooft (1986) and Tsoubelis 
(1989) and will be discussed in Section 20.4. However, the physical in- 
terpretation of such situations requires careful analysis. At this point, it 
may simply be pointed out that, in order for the matter tensor to have 
positive energy density, /' and g' must be negative on the boundaries. 



